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A COMPLEX ARITHMETIC PACKAGE 
FOR THE 3300-FORTRAN LIBRARY 



3300-FORTRAN does not do complex arithmetic. It does, 
however, provide for "type-other variables". These variables 
may have any definition that a programmer wants, provided sub- 
routines to perform the necessary operations are provided. The 
complex arithmetic package (CAP) makes use of these features of 
3300-FORTRAN. Refer to pp. C-3 and C-4 of CDC publication 
No. 60057600, Rev. B, for an explanation of the use of type- 
other variables. 

A complex number in CAP is a pair of floating-point numbers 
occupying consecutive storage locations. The first two words 
are the real part and the second two words are the imaginary 
part. All complex numbers must appear in a TYPE statement. So 
must all functions whose returned values are complex. Some 
examples are: 



TYPE COMPLEX (4) 
TYPE CO (4) 



X,Y,I,CEXP,CLOG 
J3,CC0S,CSIN,CC0SH,CSINH 



All the necessary routines to load, store, negate and mix modes 
in standard arithmetic expressions are provided except for the 
complex- to-character subroutines . The entry points in CAP are : 



R.PART 
CO.EXOR 
CO.DVOI 
CO.MUOO 
CO.MURO 
CO.SBIO 
CO.ADOR 
CO.STOI 
CO.CMOO 



I. PART 
CO.EXIO 
CO.DVOR 
CO.MUOI 
CO.SBOC 
CO.SBRO 
CO.ADIO 
CO.STOR 
CO.LNOO 



CO.EXOO 
CO.EXRO 
CO.DVIO 
CO.MUOR 
CO.SBOI 
CO.ADOO 
CO.ADRO 
CO. ST 10 
CO.LDOO 



CO.EXOI 
CO.DVOO 
CO . DVRO 
CO.MUIO 
CO . SBOR 
CO.ADOI 
CO.STOO 
CO.STRO 



ARCTAN 


COMP 


EXPI 


N0RM2 


ccos 


COSH 


IM 


XN0RM2 


CCOSH . 


CSIN 


IMAG 


RAD 


CEXP 


CSINH 


XIMAG 


RE 


CLOG 


CSQRT 


NORM 


REAL 


COM 


DIR 


XNORM 


SINH 


CONJ 









All the entries in the first group (ones with . in them), 
are standard except R.PART and I. PART. These two entries are 
actually the complex accumulator. The rest of the entries 
are functions to facilitate programming. These functions are 
described in the rest of this report. 

ARCTAN 

This is a real valued function of two real arguments . 
The results are ATAN (ARG2/ARG1) in radians, but in the 
proper quadrant. 



CCOS 



This is the complex cosine of a complex argument. 



CCOSH 



CEXP 



CLOG 



See CSINH 



This is the complex exponential of a complex argument 



This is the complex log of a complex argument. 
CLOG (X) =ALOG (SQRT (RE (X) **2+IM (X) **2) ) +COM (RAD (X) ) 



COM 



This must be declared COMPLEX. COM(ARG) leaves in the 
complex accumulator a number in the form O+iARG. ARG 
should be REAL. 



CONJ 

This gives the conjugate of a complex number. It is 
equivalent to COMP (RE (X) ,-IM (X) ) in function. 

COMP 

This must be declared COMPLEX. A function of two real 
arguments COMP (ARG1,ARG2) leave in the complex accumulator 
ARGl + iARG2. . 

COSH, SINK 

These are the cosh and sinh of a real argument. Their 
values are real. These functions are included because 
they are often useful when working with complex arithmetic. 

CSIN 

This gives the complex sine of a complex argument. 

CSINH, CCOSH 

These are the complex SINH and COSH of a complex argument. 

CSQRT 

This is the complex square root of a complex argument. 

DIR 

See RAD 

EXPI 

This is a complex function of a real argument. EXPI (THETA) 
leaves EXP(iTHETA) in the complex accumulator. 

IM, Ii4AG, XIMAG 

The three names are entry points to the same routine. This 
routine leaves in AQ the imaginary part of a complex argu- 
ment. IM and IMAG must be declared REAL. 

IMAG 

See IM 



NORM, XNORM 

These names are for computing SQRT (N0RM2 (ARG) ) . NORM 
must be declared REAL. 

N0RM2 , XN0RM2 

These names are for a routine that leaves the sum of 
the squares of the real and imaginary parts of a complex 
argument in AQ. N0RM2 must be declared REAL. 

RAD, DIR 

These are names that refer to a routine for computing the 
direction of a complex number in radians . Thus 
RAD(X)^ARCTAN(RE(X) ,IM(X) ) . See ARCTAN for a description 
of its procedure and results. 

RE, REAL 

These names are for a routine that leaves the real part 
of a complex argument in AQ. 



See 


RE 


SINK 




See 


COSH 


XIMAG 




See 


IM 


XNORM 




See 


NORM 


XN0RM2 




See 


N0RM2 



CAP is currently available under OS- 3 in FTNLIB, the 
FORTRAN library. 



